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Abstract 

In this manuscript, for approximation of solutions to equations that are 
nonlinear, a new class of two-point iterative structure that is based on a 
weight function involving two converging power series, is developed. For 
any method constructed from the developed class of methods, it requires 
three separate functions evaluation in a complete iteration cycle that is 
of order four convergence. Also, some well-known existing methods are 
typical members of the new class of methods. The numerical test on some 
concrete methods derived from the class of methods indicates that they are 
effective and competitive when employed in solving a nonlinear equation. 
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1 Introduction 


One famous iterative method (IM) that is effective when used for the ap- 
proximation of the solution s* of the nonlinear equation (NE) f(s) = 0 is the 


*Corresponding author 


Received 23 January 2022; revised 12 March 2022; accepted 2 July 2022 


Oghovese Ogbereyivwe 

Department of Mathematics and Statistics, Delta State University of Sci- 
ence and ‘Technology, Delta State, Nigeria. e-mail: ogho2015@gmail.com, 
oghovese.og bereyivweQdspz.edu.ng 


Osamede Izevbizua 
Department of Mathematics, University of Benin, Benin City, Nigeria. e-mail: 
izevbizua@uniben.edu 


Ogbereyivwe and Izevbizua 158 


Newton Iterative Method (NIM) [18], expressed as 


Wk+1 = Sk — U (Sk), (1) 
where u (sk) = tee, k = 0,1,2, ..., and |f’ (s)| 4 0. The convergence order 
(CO) of the NIM is 2 with efficiency index (EI) 1.4142. Since the emergence 
of the NIM, the plethora of techniques have been employed to modify it with 
the motivation of enhancing its CO or EI when utilized to solve an NE. One 
of such techniques that have been used is the weight function, which is a 
power series based on the form: 


w (0(s))=1+ p a (8(s))' (2) 
w (0(s))=1+ a ai (u(s) — 1)’, (3) 


where a; E€ R, 1 < m < œ and 6(s),(s) are often quotients of two real 
functions involving the evaluations of f (-) or f’ (-). 


In order to acknowledge some good literature in this direction, we note 
that, in the work of Khattri and Abbabandy [9], the NIM was adjusted and 
used as a predictor and a corrector iterative function that involves the series 
(2) to construct a one-parameter family of IM with CO four presented as 


y(Sk) = Sk — Zu(sx), 
Ski1 = Sk — U(Sk) È + De ri(v(s))*] : (4) 
n=- a, T2=(-$-304), 13 = (Ë -3a4), 


where v(s) = Fie) and agis a real parameter that is free. 


In a similar manner, the power series (2) was utilized in Chun [4] to 
suggest the two-point IM 


Wk = Sk — u(sz), 


f (we) 


Sk+1 = Sk Flsk) 


1+2 


F(we) , (E)N (5) 
Flsk) Gay | 


Ghanbari [6] used a power series like the one of (2) given as 


Af (sx)? + Bflwrk)f(sk) + eee f (we) 
Df (sx)? + Ef (we) f (se) + PF (we)? | f'(8n) 


To develop a class of three-free-parameter CO three, the IM is obtained as 


(6) 


Sk+1 = Wk | 
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1+(2+a)t,+ot?] f (we) 
L+at,+ Bt, J F(s) 


(7) 


Skil = Wk 
where t = foi) and a, 8,0 € R. It was also shown in [6] that the class of 
methods (7) has the methods presented by Ostrowski [14], Chun [4, 3], Kou, 
Li, and Wang [10], and Chun and Ham [5] as concrete members. 


In the work of Sharma and Bahl [16], the combination and power series 
of the form (2) were utilized to construct a CO four method as 


1 9f'(sk) 3f (yk)| F (sk) 
+ erty tare Pee 


(8) 


Ski1 = Sk | 


Ogbereyivwe and Ojo-Orobosa [13] employed a kind of the power series (2) 
and an additional other weight function to put forward an iterative structure 
of the form: 


1+ D (t, u) 


i=1 


G(u), (9) 


Sk+1 = Wk — 


where t = oa and ©(t,u) = ++. This iterative structure was utilized in 


1+u’ 
developing a three-parameter family of IM of CO four. 


In the separate work of Ahmad [1] and Babajee [2], the power series (3) 
was used as a weight function to develop IMs with CO four presented as 


TE sled D Sele D, (10) 


4 8 


where u(sk) = Fu) 


Also, Mahdu [11] used (3) to develop methods that utilize the same func- 
tion evaluations to achieve order four presented as 


(u(sx) = 1) 
4 


sea = se = us) G= als) [1 - +D]. a 


Motivated by the above literature that uses the power series of the kinds 
(2) or (3) as weight functions in developing some good IMs for obtaining s* 
of NE, we put forward a new class of power series-based IM that has some 
existing class of IMs that are also power series based, as concrete members. 
The method uses the quotient of two power series in the form of (3) as 
weight functions as against one power series used in Ahmad [1] and Babajee 
[2]. This technique led to the development of a class of three real and free 
parameter IM for solving NE with CO four. The manuscript is organized in 
the following sequence. Section 1 discusses some literature background, while 
Section 2 puts forward the methodology. In Section 3, the convergence test 
on the method is carried out. Section 4 presents the results of a numerical 
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experiment on the developed method. The final section of the work includes 
the conclusion and suggested areas for further research. 


2 The method 


Consider the IM (iterative method) that is based on the terminal and con- 
vergent power series put forward as: 


(: + >, ailu(sk) — n) / (: + » Bi(u(se) — | 
= m (12) 


Sk+1 = Sk — U(Sk) 


where u(s,) = a. pu (sp) = oa a; and 8; i = 1,2,3,..., are free and 


real parameters to be established. To illustrate the contribution of the IM 
(12), the case of m = 3 is considered for investigation and its convergence 
conditions. First, the following definitions are acknowledged. 


Definition 1 (Asymptotic error, Asymptotic constant, and CO). Let A; = 
sj — 8* be the IM error at the jth iteration, and suppose that the equation 
Aja = nA + O(Ag**) can be obtained from an IM by the use of the Taylor 
expansions on the functions f(-). Then Aj4+1 is called the Asymptotic error 
equation, 7 is Asymptotic constant, and p is CO. 


Definition 2 (Efficiency). If the equation Aj+1 = Af +0(AS**) is derived 
from an IM as described in Definition 1, then the value Eeff = pr (where 


T is total number of distinct functions f(-) in one IM cycle) is called the 
efficiency index of the IM. 


3 Method convergence investigation 


The investigation on the convergence of the method put forward in (12) for 
the case m = 3 is presented in this section. Firstly, the theorem that enables 
us the establishment of the method’s convergence is stated and then followed 
by its proof. 


Theorem 1. Define f: D C R > R as a function that is real valued such 
that it is differentiable for at least three times in the domain D. Suppose 
that s* C D and that |f’(-)| 4 0 in D. Then for an initial guess sọ close to 
s* in D, the sequence {8;} 59, (sj € D) of approximations, generated by the 
class of IM in (12), where m = 3, converges to s* with CO four whenever the 
free parameters a; and f; satisfy the following conditions: 8; — a, + 3 =0, 
b2 — TH (16a + 12a, — 19) = 0 and b3 — a (64a3 + 48a2 — 36a; + 54) = 0. 
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Proof. Suppose in the Taylor expansion of f (s) and f’ (s) about s* that we 
set s=s,. Then the following expression can be obtained: 


4 
f (83) = f(s") (Ay +) cn AP + 0 (43), (13) 
and 
5 
F (83) = F (s*)(1+ 50 nc, AP-*+0(A8)), — 7=0,1,2,..., (14) 


where cn = toe. n> 2. 


Using (13) and (14), the series expansions of u(s,) and ju(s,) are, respec- 
tively, 


u(sp) = Ay —c2 Aj, + (265 = 2c3) A} + (—4c3 + TC2C3 = 3c4) A% +0(4), (15) 


and 


4 8 
ulsk) =1 4 3 C2Ak (12 Sc) Az, + (36c3 — 45c2c3 + 13c4) AR 


(16) 
80 148 48 
+ ( 3 CS 3 cc3 4 77“ 5s) A% + O(A8). 
From (15), we have 
2 
Yk = Ár — zulk) 
E E AR 2 (48 —7 3c4) A% + O(A}. 
= Ae + zeA 3 (C2 c3) k T g (462 — (C23 + 3¢4) Ak + (Aj). 
(17) 
Using (17), the expansion of the function f(y,) is obtained as 
1 7 
fuk) = gfe + rhea 
Ack A(c2—c3) eN 43 
A 
( 9 3 1 a7) 7k 
4 8 2 2 
(ce (54 9 (2 ca) 9 62°83 + 5 (4c3 — 7e2c3 4 3c) A% 
a O(Ap), 
(18) 


and using (18), its derivative f’(y,) is 
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2 
f' (ye) =1+ 3024k 
1 8 4 4 
32 5°) A? 4 ( 3 (cS c3) + z203 + e) A} 


4 8 8 4 
( 3° 3 (3 a) ) C3 + 9 (204 + 3° (4c3 TC2C3 + 3c) AY 
O 


(19) 


From the expressions in (14) and (19), the Taylor expansion of the quo- 
tient £ (ue) is 


Ff’ (sx) 
Ff’ (ye) 4 
=] + -c&2 AÅ 

Tiad + 3° k 

+ [424 x A? J (36c3 45c2c3 + 13c4) A? 

273 k~ 99 2 k (20) 
80 148 32 484 
+ (F4- 3 —— eg + 3 c3 + a7 C204 5es) A% 
+ O(A}). 


Using (16), the expansion of the quotient of the two power series in (12) 
is obtained as 


3 3 
(: + », ailu(sk) — n) j (: + 3 Bilu(sk) — n) 


4 
=l-3 (a1 — 1) co Ak 


4 

Ta (4azc3 — 981c + Bi cs — 4b2c3 + ay ((9 — 481) c3 — 6c3) + 681c3) Az 
8 

T zz 88e — 368103 + 368?c3 — 867c3 + 

+ a; (4 (9 — 981 + 287 — 282) c3 + 3 (—15 + 881) c2c3 + 13c4)) A} 


1 
(1728a3c4 + 256a4c4 — 21608,c5 — 256a038,c5 + 36008?c3 + 


a 
4 (—363 + 20881) coca + 9 (32 (—3 + 281) 3 + 45c3)) At + O (AR). 
(21) 


Now, substituting (15) and (21) into the iterative structure in (12), the 
following error equation is obtained: 
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4 
Sk41 =s* + (ce + 3 (a, oa ‘i)en) A? 


+ =((-9 — 8a2 + 8a; (bı — 3) + 246, — 8 (B? — B2)) $ 
a 3 (3 + 4a, —_ 461) c3) A3 


1 
+ 57 ((4(27 + 84a2 + 16a3 — 11763 — 16023, + 848? — 168? 


+ a (117 — 846; + 166? — 1682) — 8482 + 32882 — 16f3)c3 
+ 3 (—63 — 64a + 1688; — 6487 + 8a, (881 — 21) + 6481) c2c3 
+ (81 + 104a; — 10461) c4)) AZ + O (AR) - 


(22) 

For the method (12) to attain order four convergence, the second and 

third terms of (22) must be annihilated. This can be achieved when the 
following set of equations are satisfied: 


3 
By =ar+ 1 
1 
B2 = T (16az2 + 120; — 9) (23) 
1 
Bs = Gq (Glas + 48o — 3601 +54). 


When the conditions in (23) are substituted in (22), the error equation is 
obtained as 


Aus = S” (12 = =) At +0 (AB). (24) 


Following Definition 1, the asymptotic error equation for the method (12) 
is the expression on (24), and consequently, the CO is four. 


Remark 1. By substituting the conditions (23) into (12), the class of three- 
free-parameter IM for solving NE is obtained as 
ee 14+ 07.1 ai (M(sn)—1)* 
Sk+1 = Sk — U (Sk) Essences 
B=a +3, Bq = ataa, (25) 
B — 64a3+48a2—36a01,+54 
a= 64 ; 


We emphasis here that, for any free choice of a; € R in (25), a CO four 
method with Eeff = 1.5867 can be constructed. 


Remark 2. Obtaining suitable replacement function(s) for the weight func- 
tion(s) used in [17, 13, 12, 15], and many other existing methods that were 
developed via the weight function(s) technique, requires subjecting the re- 
placement function(s) to satisfy certain conditions that involve its evaluation 
and many times differentiation at some fixed points. These procedures are 
tedious and cumbersome compared to the proposed three-parameter-based 
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class of methods (25) that requires the arbitrary substitution of values for 
the parameters. In fact, for any real values for the parameters, the method 
works. 


Remark 3. It is important to note that for m = 2 in (12), a subclass of 
method (25) is obtained as 


1+ a1 (u(se) — 1) + a2 (u(sr) — 1)? 
Sk+1 = Sk—U (Sx) 3 | (16a2+1201—9 2 j? 
1+ (a1 + 3) (u(sk)—1)+( 16 ) (u(se) — 1) 
(26) 
Some existing methods and their variants are particular members of the 
class of methods in (26). 


For instance, the method developed by Chun [4] is a concrete method of 
(26). 

Furthermore, when a; = —3 and ag = 3 in (26), the CO four methods 
by Ahmad [1] and Babajee [2] are rediscovered as 


PEET E (1 ® (lst) —1) +2 (als) nê). (27) 


Again, for ay = 3 and az = 0 in (26), the well-known Jarratt method in 
[8] is rediscovered as 


oan, n 


(u(sk) — 1) 
For some new IMs of the kind (26), the following cases are presented: 


Case 1 
Ifa, = -3 and az = 0, then the following iterative structure is discovered: 


1- $ (Msk) — 1) ) 
1 — $ (u(s) — 1)" 


bolo} ew 


Sk+1 = Sk — U (sp) (; 


Mı: Sk+1 = Sk — U (Sp) ( 


Case 2 
The substitution of a; = + and az = 4 in (26) yields the method 


M» : unsa un (HEE D tte i. 


1+ 34 (alsk) — 1) + £ (uls) — 1)? 


Case 3 
Consider a, = 4 and a2 = 0. Then a new method is obtained as 


1 
M3:  Sk41 = Sk — uU (sp) ( 
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Remark 4. In order to illustrate some new methods of kind (25), the fol- 
lowing IMs are constructed for some real values of a;: 

Case 4 
When a, = a2 = a3 = 0, then the method My, is constructed as 


Sets = 8h — UES) k e-d 


Case 5 
Consider the case ay = ag = a3 = 3. Then the IM M; is obtained as 


= sp — U (Sk) 1+3 ((u(sn) 1) + (u(sk)— 1) + (44(se) 1)°) 
ee SNES (ulsx) — 1) + 8 (usu) — 1? + 2 (use) — 


4 Numerical test 


In this section, the methods derived from the developed class of methods are 
subjected to numerical tests in order to verify their effectiveness when used 
to obtain a solution to NE. The numerical results of the new methods are 
compared with the results of some existing methods that are of CO four. 
These include the method developed by Sharma and Bahl [16] (presented in 
(8)) denoted by SBM, method developed by Ogbereyivwe and Ojo-Orobosa 
in [13] (presented in (9)) with the parameters 8 = 2,7 = ø = 0, denoted 
by OOM, the method of Hafiz and Khirallah [7], denoted as HKM, the bi- 
parameter method in Shams et al. [15] (given in (28)) with 8 = 1 and a = 2 
and denoted as SEM, the methods in Ghanbari [6] (presented in (7)) obtained 
by taking a = —2, 8 = 1,0 = 4 and a = —2,8 = 1,0 = —1, and denoted as 
GM1 and GM2, respectively. The metric used for comparison includes the 
computational CO peoc is given by 


logio If (se+a)I/ [F (se) | 
logio |f (s)| / |f (8x1) 


and |f (s,)|. All computation programs for the methods were written and 
executed in the MAPLE 2017 software environment using Intel Celeron(R) 
computer with a 2GB RAM processor. The error tolerance |f (s)| < 2007200 
was used in terminating program execution. The computational outputs are 
maintained at 2000 digits precision in all computations. The NE f;(s) 
0 used for the test and their respective solutions s* presented in 15 digits 
precision are as given in Examples 1—5. For comparison, computation results 
are presented in Tables 1 and 2. 
Consider the following examples: 


Peoc = 


Example 1. [6] fi(s) = s?—e*—-3s+2 = 0, s* = 0.257530285439860.... 
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Example 2. [6] fo(s) = sin?(s)—s?+1 = 0, s* = 1.404491648215341.... 
Example 3. [13] f3(s) = arctan(s) = 0, s*=0. 

Example 4. [6] fu(s) =(s—1)*>-2=0,  s* = 2.259921049894873.... 


Example 5. [13] fa(s) = f+ 2 logy (4 + 251) = 0, where $ = 1074 
and R = 10°, s* = 0.004129319146099.... 


Table 1: Methods results comparison for Examples 1- 2 


Methods so JED VGD VGI VCI Fs) Pa 
SBM 5.3e—4 7.le—18 2.3e—73 2.7e — 295 — 4.0 
HKM T.3e—3 3.4e—16 1.6e—66 7.le— 268 - 4.0 
SEM 5.9e—3 1.7e—13 1.le—55 1.7e— 224 - 4.0 
OOM 1.7e—2 2.9e—11 2.5e—46 1.4e— 186 - 4.0 
GM1 2.8&e—4 3.4e—19 7.5e—79 1.7e— 317 - 4.0 
GM2 -0.3 3.0e—3 1.le—14 1.6e—60 8.4e — 244 - 4.0 

M1 1.6e—4 4.8e—?20 4.8e—82 4.3e — 330 - 4.0 
M2 2.3e—3 3.9e—15 3.le—62 1.2e — 250 - 4.0 
M3 T.3e—5 2.3e—21 2.5e—87 3.le— 351 - 4.0 
M4 2.le—3 2.6e—15 5.le—63 1.5e — 253 - 4.0 
M5 2.4e—3 4.0e—15 3.3e—62 1.6e — 250 - 4.0 
SBM 5.3e—2 45e€-7 2.8e—27 4.0e—108 1.8e— 431 4.0 
HKM 3.7Te—2 6.2e—8 5.0e—31 2.le— 123 6.3e— 493 4.0 
SEM 1.2e—2 1.5e—9 4.1e—37 2.2e—147 1.7s—588 4.0 
OOM 1.0e—2 1.6e-5 1.7e —20 2.2e—80 3.6e—320 4.0 
GM1 6.le—2 1.0e—6 1.le—25 1.le— 101 1.4e— 405 4.0 
GM2 2 1.4e—2 2.8e—9 3.9e—36 1.6e— 143 4.7e—573 4.0 
M1 5.6e—2 7.5e—7 3.0e—26 6.7e—104 1.8e—414 4.0 
M2 6.le—2 1.2e—6 1.9e—25 1.3e— 100 2.9e—401 4.0 
M3 5.9e—2 9.3e—7 68e-—26 1.9e— 102 1.le—408 4.0 
M4 2.7e—2 2.le—9 9.8e—38 5.le—151 3.6e—604 4.0 
M5 1.5e—2 2.4e—10 1.9e—41 6.8e—166 1.le—663 4.0 


From the numerical results in Tables 1 and 2, we observe that all the par- 
ticular methods constructed from the developed class of methods (25) solved 
the test problems very well. Furthermore, the new methods’ computational 
order of convergence displayed in the last column of Tables 1 and 2 agree 
with the theoretical CO established in Section 3. The methods computation 
results outputs are also highly competitive with the methods compared. 


5 Conclusion 


The new three-parameter class of IM put forward in this manuscript was 
shown to have a capacity of approximating the solution of NE with CO 
four requiring three evaluations of distinct functions in an iteration cycle. 
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Table 2: Methods results comparison for Examples 3- 5 


Methods so  |£F(s)|  |£F(s2)| |f (ss) If (s4)| [F(s5)| Peoe 
SBM 6.8e—3 2.3e-— 12 1.le — 59 2.2e — 296 - 4.0 
HKM 4.7Te—3 3.7e—13 1.le — 63 2.3e — 316 - 4.0 
SEM 4.7Te—3 5.0e-—13 6.8e — 62 3.2e — 312 - 4.0 
OOM 6.le-2 15e-7 1.6e — 35 2.6e — 175 - 4.0 
GM1 9.3e—-3 1.6e— 11 2.2e — 55 1.le — 274 - 4.0 
GM2 —0.5 13e-3 7.3e-— 16 4.6e — 77 4.8e — 383 - 4.0 

M1 9.2e-3 1.le—11 2.4e — 56 1.2e — 279 4.0 
M2 2.9e—3 3.le—14 4.8e — 69 4.4e — 343 - 4.0 
M3 3.le—3 4.4e-— 14 2.7e — 68 2.2e — 339 - 4.0 
M4 1.5e—-3 14e-15 9.6e — 76 1.3e — 376 - 4.0 
M5 3.le—3 4.5e-— 14 3.2e — 68 5.6e — 339 - 4.0 
SBM 2.5e—1 2.5e-—5 3.8e — 21 2.0e— 84 1.3e—337 4.0 
HKM 1.5e—1 1.8e-—6 4.8e — 26 2.2e— 104 1.0e— 417 4.0 
SEM 8.2e-1 7.2e-2 4.0e — 7 3.0e— 28 1.0e—112 4.0 
OOM 4.6e—1 7.5e—4 1.0e — 14 4.0e— 58 8.5e—232 4.0 
GM1 2.9e—-1 57e-—5 1.3e — 19 3.3e— 78 1.4e— 312 4.0 
GM2 3 3.2e—-2 1.le-—8 1.5e — 34 5.9e — 138 1.3e— 551 4.0 
M1 2.7e—1 42e-—5 3.8e — 20 2.6e— 80 5.8e—321 4.0 
M2 3.5e—-2 2.3e-—9 4.4e — 38 5.8e — 153 1.7e—612 4.0 
M3 3.4e—2 6.6e— 10 1.0e — 40 5.2e — 164 3.8e— 657 4.0 
M4 9.6e—-2 94e-8 1.2e — 31 3.5e— 127 2.2e— 509 4.0 
M5 1.2e—-2 3.3e—11 1.9e — 45 2.le— 182 3.2e—730 4.0 
SBM 3.7e—1 1.4e-5 2.6e — 23 3.7e—94 1.5e—377 4.0 
HKM 1.6e—3 8.5e—16 7.6e — 65 4.3e — 261 - 4.0 
SEM Failed to converge 
OOM 6.le—-1 34e-—4 2.8e — 17 1.4e— 69 8.7e-—279 4.0 
GM1 4.5e—1 3.7e-—5 2.0e — 21 1.6e— 86 6.8e—347 4.0 
GM2 0.002 2.9e-1 1.5e-5 8.0e — 23 7.3e—92 5.2e— 368 4.0 
M1 4.le—1 29e-5 7.8e — 22 4.3e— 88 4.0e—353 4.0 
M2 1.0509 8.le— 4 2.4e — 16 1.9e— 66 7.4e—267 4.0 
M3 5.le—1 3.3e-5 4.9e — 22 2.6e—89 2.0e—358 4.0 
M4 1.2e—1 14e-7 1.9e — 31 8.le— 127 24e—508 4.0 
M5 7.3e—1 17e-4 5.0e — 19 3.6e—77 9.7e—310 4.0 


Several famous existing methods, such as Chun [4], Ahmad [1], Babajee [2], 
Jarrat [8], and others, are typical members of the class of developed methods. 
For future research, the implementation of the method developed herein in 
Banach spaces may be considered. 
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